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ABSTRACT 



In this talk we will discuss three different issues. First of all, there exist several 
proposals how to solve cosmological problems by adiabatic expansion of the Universe 
of a specific type, without any use of inflation. We explain why these models do not 
solve the flatness/entropy problem. On the other hand there exist some claims that 
inflation also does not solve the flatness problem. We show that these claims are 
correct to a certain extent for old and new inflation, but are incorrect for simplest 
versions of chaotic inflation. 

The second issue to be considered is the classification of inflationary models. At 
present all kinds of inflationary models are compared to each other in order to un- 
derstand which of them better describes observational data. Unfortunately, this often 
leads to terminological misunderstandings and even to misleading conclusions. We 
discuss three coexisting types of classification of inflationary models and explain their 
relation to each other. We also briefly describe the hybrid inflation model, which 
combines several interesting features of different inflationary models. 

Finally we show that it is possible to describe inflationary universe by a stationary 
(i.e. time- independent) probability distribution to find a domain of a given size con- 
taining matter with given properties. This represents a strong deviation of inflationary 
cosmology from the standard Big Bang paradigm. 



1 Based on the talks given at the 5th Canadian Conference on General Relativity, Waterloo, at the Yamada 
Conference on Cosmology, Tokyo, and at the Summer School on Cosmology in Aspen, 1993 
2 On leave from: Lebedev Physical Institute, Moscow. E-mail: linde@physics.stanford.edu 



1 Introduction 



After more than ten years of its existence, the inflationary universe scenario is gradually becoming 
a standard cosmological paradigm. Inflation solved many different problem of the standard 
cosmology, including the horizon, flatness, homogeneity and isotropy problems. For a detailed 
discussion of this theory see e.g. [|TJ. However, it would be incorrect to stop looking for other 
solutions of these problems. Moreover, one should remember that not every inflationary model 
solves these problems well enough. 

Recently Freese and Levin M proposed a model which they called MAD. According to [0, 
strong variation of the gravitational constant in the very early universe may solve the horizon 
problem without any need for inflation and for any nonadiabatic processes such as reheating of 
the Universe after inflation. Their model has some interesting features, but it is not clear at all 
whether this model can solve the homogeneity and isotropy problem, and it appears that MAD 
scenario cannot solve the flatness problem. 

A detailed discussion of the last issue was contained in a paper by Hu, Turner and Weinberg 
0. However, this discussion was concentrated on the model of ref. |2j and on its simplest 
generalizations. Meanwhile we all remember that the version of the inflationary universe scenario 
suggested by Guth Q also did not work, and a year later Guth and Weinberg have almost 
completely proved that a consistent inflationary theory is impossible 0. What if a year from 
now somebody will suggest a new or extended MAD scenario? 

Unfortunately (or fortunately, depending on the point of view), there are some reasons to 
believe that this is impossible: all models assuming adiabatic (or almost adiabatic) expansion of 
the universe, including the MAD scenario |J, have a serious problem "by definition". Indeed, 
the total entropy of the observable part of the Universe is greater than S ~ 10 88 . If expansion 
of the Universe is adiabatic, one is forced to assume that the Universe from the very beginning 
was huge: it contained more than 10 88 particles. This is a manifestation of the entropy problem, 
which is almost equivalent to the flatness problem, but is more general [|TJ. 

The best way to understand that this problem is more general than the flatness problem is 
to remember the proposal to use Kaluza-Klein compactification as an alternative to inflation, 
see e.g. ||. The main idea is to begin with a Universe of a very large number of dimensions, 
say, D = 100. Let us assume that the Universe initially was closed, and its initial size Iq was 
just few times greater than the Planck length l p ~ M~ x . However, this is enough to incorporate 
enormously large amount of particles into such a universe, because of its huge volume. At the 
Planck temperature T ~ M p such a universe contained ~ (ZoT) 99 particles, which is a very large 
number if the universe originally was few times greater than the Planck length. As a result of 
compactification, all these particles should fit into a Universe with only three spatial dimensions, 
which leads to its extremely rapid expansion. 

Thus, in this scenario the Universe is initially very curved. Then, due to compactification 
of almost all its dimensions, the Universe becomes three-dimensionally flat without entering a 
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stage of vacuum energy dominance and subsequent reheating. One may argue that this solves 
the flatness problem without any need of inflation. However, one should still assume that the 
Universe contained at least 10 88 particles from the very beginning, which is a manifestation of 
the entropy problem. Until this problem is unresolved, all models which do not assume strong 
violation of adiabaticity cannot compete with inflation. 

But is it correct that inflation solves this problem? The answer depends on the choice of 
inflationary scenario. Let us consider for example a closed universe in the new inflationary 
universe scenario, which occurs when the energy density of ultrarelativistic particles becomes 
smaller than the effective potential V(0). One can easily show that this happens at the time 




2vrU(0) ' (1) 

Meanwhile, as it is shown in Jl|, a closed universe filled by the ultrarelativistic gas with a total 
entropy S begins collapsing after the time 

1/6 

S 2/3 M- 1 . (2) 
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If this happens before the Universe becomes vacuum energy dominated, then the Universe recol- 
lapses and no inflation occurs. Using eqs. (|IJ) and (^j) one can easily check that inflation in this 
model is possible only if 

For a typical potential which appears in grand unified theories V(0) ~ 10~ 15 Mp. This requires 
S > 10 11 . This is an improvement by at least 77 orders of magnitude as compared with any 
adiabatic expansion scenario, including the MAD model || or the Kaluza-Klein scenario ||. 
However, 10 11 is still a very big number. At the Planck time the total mass of the Universe 
containing 10 11 particles with the Planck energy M p ~ 10 19 GeV ~ 10~ 5 g would be about ~ 10 6 

g ~ 10 M Schwarzenegger^ which is a lot! 

Fortunately, eq. (|3|) suggests an easy resolution of this problem. If inflation may begin at a 
density not much lower than the Planck density, the entropy problem disappears. This is precisely 
the case in the simplest models of chaotic inflation based on the potentials ~ (f) n or ~ e°"^. In 
such theories inflation may begin in a state with a density close to the Planck density and with 
vanishing entropy. This solves both the flatness and the entropy problem M. 



3 We use the system of units M Schwarzenegger = 10 5 g, which differs by ten orders of magnitude from the Planck 
mass M p — 10~ 5 g. 
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2 Classification of inflationary models and hybrid Infla- 
tion 



After the discovery of anisotropy of the microwave background radiation by COBE, there were 
many papers trying to compare observational consequences of various inflationary models. How- 
ever, some of the authors of these papers used their own definitions of new inflation, chaotic 
inflation, etc. This gave rise to certain purely terminological misunderstandings. For example, 
some authors believe that if chaotic inflation is right then extended inflation is wrong (and vice 
versa), that the difference between chaotic inflation and the new inflation is that the new inflation 
begins at small 0, whereas chaotic inflation begins at large <fi. These and some other misconcep- 
tions appear due to the existence of several different classifications of inflationary models. 

The first classification deals with the initial conditions for inflation. The old and new infla- 
tionary models were based on the assumption that the Universe from the very beginning was in 
a state of thermal equilibrium at an extremely high temperature, and that the inflaton field 
was in a state corresponding to the minimum of its temperature-dependent effective potential 
V((f>) j|, 0| . The main idea of the chaotic inflation scenario was to study all possible initial condi- 
tions in the Universe, including those which describe the Universe outside of the state of thermal 
equilibrium, and the scalar field outside of the minimum of U(0) ||. This scenario includes 
the possibility of inflation from the state with = in a thermal equilibrium, like in the new 
inflationary Universe scenario, but it contains many other possibilities as well. Therefore it can 
be realized in a much greater variety of models than the new inflationary Universe scenario. In 
fact, at present the idea of thermal beginning is almost completely abandoned, and all realistic 
models of inflation from the point of view of the first classification are of the chaotic inflation 

type !■ 

The second classification describes various regimes which are possible during inflation: quasi- 
exponential inflation, power law inflation, etc. This classification is absolutely independent of the 
issue of initial conditions. Therefore it does not make any sense to compare, say, power law infla- 
tion and chaotic inflation, and to oppose them to each other. For example, in |5| it was pointed 
out that the chaotic inflation scenario, as distinct from the new inflationary Universe scenario, 
can be realized in the theories with the effective potential e°"^ for a <C Vl67r. Meanwhile, in 
|T0|| it was shown that this inflation is power law. Thus, the inflationary Universe scenario in the 



theory e a< ^ describes chaotic power law inflation. 

Finally, the third classification is related to the way inflation ends. There are two possibilities 
extensively discussed in the literature: slow rollover versus the first-order phase transition. The 
models of the first class describe slow rolling of the inflaton field 0, which gradually becomes 
faster and faster. A particular model of this type is chaotic inflation in the theories 0™. The 
models of the second class should contain at least two scalar fields, and a. They describe a 
strongly first-order phase transition with bubble production which is triggered by the slow rolling 



of the field 0. One of the popular models of this type is the extended inflation scenario |TT| , which 



is a combination of the Brans-Dicke theory and the old inflationary scenario. There exist other 
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versions of the first-order scenario with two scalar fields, which do not require any modifications 
of the Einstein theory of gravity, see e.g. [12]. 



In the beginning it was assumed that the bubbles formed during the first-order phase transi- 
tion could be a useful ingredient of the theory of the large scale structure formation. However, 
later it was realized that one should make considerable modifications of the original models in 
order to avoid disastrous consequences of the bubble production. According to the most recent 
modification [|l3j, the bubble formation happens only after the end of inflation. In this case, the 
end of inflation occurs as in the standard slow-rollover scenario. Therefore it would be interesting 
to find out other possible ways in which inflation may end in the models with several different 
scalar fields. More generally, one may try to find out other qualitatively new inflationary regimes 
which may appear due to a combined evolution of several scalar fields. 

Of course, one should not invent excessively complicated models without demonstrated need. 
However, sometimes qualitatively different inflationary regimes appear after minor modifications 
of the basic inflationary models, or after making their hybrids. For example, in |L4|, [15|] we 
proposed a very simple model of two interacting scalar fields where inflation may end by a rapid 
rolling of the field a ('waterfall') triggered by the slow rolling of the field <fi. This regime differs 
both from the slow-rollover and the first-order inflation. By changing parameters of this model 
one can continuously interpolate between these two regimes. Therefore some hybrid models of 
such type may share the best features of the slow-rollover and the first-order models. 



The effective potential of the hybrid inflation model of ref. [14, 15 1 is given by 



V(a,<f>) 



4A 



(M - Xa 2 ) 2 + 



2^2 



2 2 

ra 



(4) 



Theories of this type were considered in JT6[-|[T8|. The main difference between the models of 



refs. [|TI|-[|]1| and our model is a specific choice of parameters, which allows the existence of the 
waterfall regime mentioned above. There is also another important difference: we will assume 
that the field a in this model is the Higgs field, which remains as a physical degree of freedom 
after the Higgs effect in an underlying gauge theory with spontaneous symmetry breaking. This 
field acquires only positive values, which removes the possibility of domain wall formation in this 
theory. Usually it is rather dangerous to take the inflaton field interacting with gauge fields, since 
its effective coupling constant A may acquire large radiative corrections ~ e 4 , where e is the gauge 
coupling constant. In our case this problem does not appear since density perturbations in our 
model remain small at rather large A, see below. 



The effective mass squared of the field a is equal to —M 2 + g 2 <p 2 . Therefore for <fi > <p c = M/g 
the only minimum of the effective potential V(a, <fi) is at a = 0. The curvature of the effective 
potential in the cr-direction is much greater than in the ^-direction. Thus we expect that at the 
first stages of expansion of the Universe the field a rolled down to o = 0, whereas the field 
could remain large for a much longer time. For this reason we will consider the stage of inflation 
at large 0, with o = 0. 
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At the moment when the inflaton field becomes smaller than C = M/g, the phase transition 
with the symmetry breaking occurs. If m 2 2 = m 2 M 2 /g 2 < M 4 /A, the Hubble constant at the 
time of the phase transition is given by 

Thus we will assume that M 2 3> ^jh. We will assume also that m 2 <C H 2 , which gives 

M 2 > mM p \ — . (6) 
V 2tc 



One can easily verify, that, under this condition, the Universe at > C undergoes a stage of 
inflation. In fact, inflation in this model occurs even if m 2 is somewhat greater than H 2 . Note 
that inflation at its last stages is driven not by the energy density of the inflaton field but by 
the vacuum energy density V(0, 0) = as in the new inflationary Universe scenario. This was 
the reason why we called this model 'hybrid inflation' in . 



Let us study the behavior of the fields and a after the time At = H~ x = yf^;j^f from the 
moment t c when the field becomes equal to C . The equation of motion of the field during 
inflation is 3H<p = m 2 (f). Therefore during the time interval At = H^ 1 the field decreases from 

C by A0 — ^jgr — 2ngM% ' ^ ne a bsolute value of the negative effective mass squared —M 2 +g 2 (p 2 
of the field a at that time becomes equal to 

9/ s \m 2 M 2 

M <*> = • (7) 

The value of M 2 (0) is much greater than H 2 for M 3 A mM 2 . In this case the field a within 
the time At ~ H^ 1 rolls down to its minimum at <r(0) = M(0)/vA, rapidly oscillates near it and 
loses its energy due to the expansion of the Universe. However, the field cannot simply relax near 
this minimum, since the effective potential V(0, a) at <r(0) has a nonvanishing partial derivative 

dV 2 g 2 <PM 2 (<p) 

50 =m0+ ^^' (8) 

One can easily check that the motion in this direction becomes very fast and the field rolls to the 
minimum of its effective potential within the time much smaller than if -1 if M 3 <^ \/~Xg mM 2 . 
Thus, under the specified conditions inflation ends in this theory almost instantaneously, as soon 
as the field reaches its critical value C = M/g. 

The amplitude of adiabatic density perturbations produced in this theory can be estimated 
by standard methods M and is given by 



dp _ WVQ^ V 3 / 2 16 v / 6^(^ + 



rn 2 <f>^ 



p 5 M 3 ^ 5M 3 m 2 



3/2 

-. (9) 
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In the case m 2 <C H 2 the scalar field does not change substantially during the last 60 e-foldings 
(i.e. during the interval At ~ 60i/ _1 ). In this case the amplitude of density perturbations 
practically does not depend on scale, and is given by 



5 P 2V^gM 5 



The definition of ^ used in Q corresponds to COBE data for ^ ~ 5 - 10 5 . Dividing it by (|6]) with 
an account taken of (fl0|) gives M 3 <C 5- 10~ 5 \g~ 1 mM 2 . This means that the 'waterfall conditions' 



M 3 <C A mM 2 and M 3 <C V\g mM 2 automatically follow from the conditions m 2 <C H 2 and 
^ ~ 5 • 10 -5 , unless the coupling constants A and g are extremely small. Therefore the waterfall 
regime is realized in this model for a wide variety of values of parameters m, M, A and g which 
lead to density perturbations ~ 5 ■ 10~ 5 . 

To give a particular example, let us take g 2 ~ A ~ 10 _1 , m ~ 10 2 GeV (electroweak scale). 

In this case all conditions mentioned above are satisfied and ^ ~ 5 • 10~ 5 for M ~ 1.3 • 10 11 GeV. 

p 

In particular, we have verified, by solving equations of motion for the fields and a numerically, 
that inflation in this model ends up within the time At <C H^ 1 after the field reaches its critical 
value C = M/g. The value of the Hubble parameter at the end of inflation is given by H ~ 7- 10 3 
GeV. The smallness of the Hubble constant at the end of inflation makes it possible, in particular, 
to have a consistent scenario for axions in inflationary cosmology even if the axion mass is much 
smaller than 10~ 5 eV This model has some other distinctive features. For example, the 



spectrum of perturbations generated in this model may look as a power-law spectrum rapidly 
decreasing at large wavelength I [ffSl ITS . 



Indeed, at the last stages of inflation (for ^ ^> the field behaves as 

/ m 2 (t — t c )\ , , 

= C ■ exp ( ^-^) , (11) 

whereas the scale factor of the Universe grows exponentially, a ~ e . This leads to the following 
relation between the wavelength of perturbations / and the value of the scalar field at the 

n yn 2 /3H 2 

moment when these perturbations were generated: ~ C ( f-j .In this case 

Sp 2VfagM 5 (lY^ 



P 5\V\M^m 2 



which corresponds to the spectrum index n = 1 + = 1 + A " A ^f p 10 -Q Note that this spectrum 



index is greater than 1, which is a very unusual feature. For the values of m, M, A and g considered 
above, the deviation of n from 1 is vanishingly small (which is also very unusual). However, let 
us take, for example, X = g = 1, M = 10 15 GeV (grand unification scale), and m = 5 • 10 10 GeV. 
In this case the amplitude of perturbations at the end of inflation (0 = C ) is equal to 4 • 10~ 4 , 
n ~ 1.1, and the amplitude of the density perturbations drops to the desirable level — ~ 5 ■ 10~ 5 

4 This result recently was also obtained in p7|. 
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on the galaxy scale (l g ~ l c ■ e 50 ). One may easily obtain models with even much larger n, but 
this may be undesirable, since it may lead to formation of many small primordial black holes [20 . 



Note, that the decrease of -j at large / is not unlimited. At > -^j- the spectrum begins 
growing again. Thus, the spectrum has a minimum on a certain scale, corresponding to the 
minimum of expression (P). This complicated shape of the spectrum appears in a very natural 
way, without any need to design artificially bent potentials. 

As we have seen, coupling constants in our model can be reasonably large, and the range of 
possible values of masses m and M is extremely wide. Thus, our model is very versatile. One 
should make sure, however, that the small effective mass of the scalar field does not acquire 
large radiative corrections near <ft = <p c . Hopefully this can be done in supersymmetric theories 
with flat directions of the effective potential. 



One can suggest many interesting generalizations of our model. For example, instead of the 
term ? in (f|) one can use the term In this case one may have two disconnected stages of 

inflation. The first stage occurs at large <ft, as in the simplest version of chaotic inflation scenario. 
This stage ends at (f> < M p /3, if M 2 <C X^M 2 . Then the field rapidly rolls down and oscillates 

until the amplitude of its oscillations becomes smaller than ~ X M ^ . At this moment the fre- 



quency of oscillations ~ v becomes smaller that the Hubble constant, and the second stage 
of inflation begins. This stage of inflation ends with the waterfall at <p c = M/g. As was shown 
in [plU , in the models with two stages of inflation with a break between them the spectrum of 
density perturbations may have a very reach and non-trivial structure. 



3 Stationary Universe 



The first models of inflation were based on the standard assumption of the Big Bang theory that 
the Universe was created at a single moment of time in a state with the Planck density, and that 
it was hot and large (much larger than the Planck scale M~ l ) from the very beginning. The 
success of inflation in solving internal problems of the Big Bang theory apparently removed the 
last doubts concerning the Big Bang cosmology. It remained almost unnoticed that during the 
last ten years the inflationary theory has broken the umbilical cord connecting it with the old Big 
Bang theory, and acquired an independent life of its own. For the practical purposes of description 
of the observable part of our Universe one may still speak about the Big Bang. However, if one 
tries to understand the beginning of the Universe, or its end, or its global structure, then some 
of the notions of the Big Bang theory become inadequate. 

For example, already in the first version of the chaotic inflation scenario || there was no 
need to assume that the whole Universe appeared from nothing at a single moment of time 
associated with the Big Bang, that the Universe was hot from the very beginning and that the 
inflaton scalar field <fi which drives inflation originally occupied the minimum of its potential 
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energy. Later it was found that if the Universe contains at least one inflationary domain of a 
size of horizon ('/i- region') with a sufficiently large and homogeneous scalar field 0, then this 
domain will permanently produce new /i-regions of a similar type. In other words, instead of 
a single Big Bang producing a one-bubble Universe, we are speaking now about inflationary 
bubbles producing new bubbles, producing new bubbles, ad infinitum. In this sense, inflation is 
not a short intermediate stage of duration ~ 10~ 35 seconds, but a self-regenerating process, which 
occurs in some parts of the Universe even now, and which will continue without end. The most 
striking realization of this scenario occurs in the context of chaotic inflation | 22|| , but the basic 
features of this scenario remain valid in old inflation 
inflation as well 



231, new inflation [24, 251 and extended 
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Thus, recent developments of inflationary theory have considerably modified our cosmological 
paradigm |IJ. Now we must learn how to formulate physical questions in the new context. For 
example, in a homogeneous part of the Universe there is a simple relation between the density of 
matter and time. However, on a very large scale the Universe becomes extremely inhomogeneous. 
Its density, at the same 'cosmic time', varies anywhere from zero to the Planck density. Therefore 
the question about the density of the Universe at the time 10 10 years may not have any definite 
answer. Instead of addressing such questions we should study the distribution of probability of 
finding a part of the Universe with given properties, and find possible correlations between these 
properties. 



It is extremely complicated to describe an inhomogeneous Universe and to find the corre- 
sponding probability distribution. Fortunately, there exists a particular kind of stationarity of 
the process of the Universe self-reproduction which makes things more regular. Due to the no- 
hair theorem for de Sitter space, the process of production of new inflationary domains occurs 
independently of any processes outside the horizon. This process depends only on the values of 
the fields inside each /i-region of radius if -1 . Each time a new inflationary /i-region is created 
during the Universe expansion, the physical processes inside this region will depend only on the 
properties of the fields inside it, but not on the 'cosmic time' at which it was created. 



In addition to this most profound stationarity, there may also exist some simple stationary 
probability distributions which may allow us to say, for example, what the probability is of finding 
a given field at a given point. To examine this possibility one should consider the probability 
distribution P c (4>, t\x), which describes the probability of finding the field at a given point at a 



time t, under the condition that at the time t = the field at this point was equal to x HI) 26 



The same function may also describe the probability that the scalar field which at time t was 
equal to 0, at some earlier time t = was equal to x- 



The probability distribution P c has been studied by many authors, see e.g. P2| , p5|-|f27f. 
Our investigation of this question has shown that in all realistic inflationary models the proba- 
bility distribution P c (0, t\x) is not stationary [ESJ. The reason is very simple. The probability 



distribution P c is in fact the probability distribution per unit volume in comoving coordinates 
(hence the index c in P c ), which do not change during expansion of the Universe. By considering 
this probability distribution we neglect the main source of the self-reproduction of inflationary 
domains, which is the exponential growth of their volume. Therefore, in addition to P c , we in- 
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troduced the probability distribution P p (<f>,t\x), which describes the probability to find a given 
field configuration in a unit physical volume [22|. In the present paper (see also |28| for a more 



detailed presentation) we will show that under certain conditions the stationary probability distri- 
bution P p (4>, t\x) does exist, and a typical relaxation time during which the distribution P p ((f>, t\x) 
approaches the stationary regime is extremely small. 



First of all we should remember some details of stochastic approach to inflation. Let us con- 
sider the simplest model of chaotic inflation based on the theory of a scalar field minimally 
coupled to gravity, with the effective potential V(4>). If the classical field <fi is sufficiently homo- 
geneous in some domain of the Universe, then its behavior inside this domain is governed by the 
equations 

<j> + 3Hj>= -dV/d(f) , (13) 

Here H = a/a, a(t) is the scale factor of the Universe, k = +1, —1, or for a closed, open or flat 
Universe, respectively. M p is the Planck mass, which we will put equal to one in the rest of the 
paper. 

Investigation of these equations has shown that for many potentials V(<f)) (e.g., in all power- 
law V((p) ~ cf) n and exponential V(<f)) ~ e°"^ potentials) there exists an intermediate asymptotic 
regime of slow rolling of the field (f> and quasi-exponential expansion (inflation) of the Universe 
U. At this stage the Hubble parameter is H((f)) = y / 87rV A (0)/3. In the theories V((f)) ~ <p n 
inflation ends at <p e ~ n/12. In the theory with V(<f)) ~ e 01 ® inflation ends only if we bend the 
potential at some point e ; for definiteness we will take <fi e = in this theory. 

Inflation stretches all initial inhomogeneities. Therefore, if the evolution of the Universe were 
governed solely by classical equations of motion, we would end up with an extremely smooth 
Universe with no primordial fluctuations to initiate the growth of galaxies. Fortunately, new 
density perturbations are generated during inflation due to quantum effects. The wavelengths of 
all vacuum fluctuations of the scalar field grow exponentially in the expanding Universe. When 
the wavelength of any particular fluctuation becomes greater than if -1 , this fluctuation stops 
oscillating, and its amplitude freezes at some nonzero value 5<p(x) because of the large friction 
term 3H<j) in the equation of motion of the field (p. The amplitude of this fluctuation then remains 
almost unchanged for a very long time, whereas its wavelength grows exponentially. Therefore, 
the appearance of such a frozen fluctuation is equivalent to the appearance of a classical field 
5(j)(x) that does not vanish after averaging over macroscopic intervals of space and time. 

Because the vacuum contains fluctuations of all wavelengths, inflation leads to the creation 
of more and more perturbations of the classical field with wavelengths greater than H~ l . The 
average amplitude of such perturbations generated during a time interval H~ l (in which the 
Universe expands by a factor of e) is given by 

|^(ar)|«^. (15) 

Z7T 
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The phases of each wave are random. Therefore, the sum of all waves at a given point fluctuates 
and experiences Brownian jumps in all directions in the space of fields. 



The standard way to describe the stochastic behavior of the inflaton field during the slow- 
rolling stage is to coarse-grain it over /i-regions and consider the effective equation of motion of 
the long- wavelength field [ 2t| : 



dt 3H(<f>) 2tt ^ lJ ' 



(16) 



Here H = y8irV/3, £(£) is the effective white noise generated by quantum fluctuations, which 
leads to the Brownian motion of the classical field 6. 



This Langevin equation leads to two stochastic equations for the probability distribution 
P c ((p,t\x)- The first one is called the backward Kolmogorov equation, 



dt ~~ 2 2vr dx\ 2vr d\ 



V'ix) d 



P c {<f>,t\x) 



(17) 



In this equation one considers the value of the field at the time t as a constant, and finds the 
time dependence of the probability that this value was reached during the time t as a result of 
diffusion of the scalar field from different possible initial values x = 0(0)- 

The second equation is the adjoint to the first one; it is called the forward Kolmogorov 
equation, or the Fokker-Planck equation, 

dP c {(f),t\x) 1 d /# 3/2 (0) d /# 3/2 (0) 



dt 



2tt 



This equation was derived in [26|, see also [25 



One may try to find a stationary solution of equations (|17D , (fQj), assuming that 9Pc ^ t ^ = 0. 
The simplest stationary solution (subexponential factors being omitted) would be 



PciMx) ~ ex P 



• exp 



(19) 



V sv (x)J ^ 

This function is extremely interesting. Indeed, the first term in (|l9|) is equal to the square of the 



Hartle-Hawking wave function of the Universe [29|, whereas the second one gives the square of 



the tunneling wave function 30 1 ! 



At first glance, this result gives a direct confirmation and a simple physical interpretation 
of both the Hartle-Hawking wave function of the Universe and the tunneling wave function. 
However, in all realistic cosmological theories, in which V(4>) = at its minimum, the Hartle- 
Hawking distribution exp (gp 3 ^) is not normalizable. The source of this difficulty can be easily 
understood: any stationary distribution may exist only due to compensation of the classical flow 
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of the field downwards to the minimum of V(0) by the diffusion motion upwards. However, 
diffusion of the field discussed above exists only during inflation. Thus, there is no diffusion 
motion upwards from the region < e . Therefore all solutions of equation flUf ) with the proper 
boundary conditions at = e (i.e. at the end of inflation) are non-stationary (decaying) [ p2]| . 



The situation with the probability distribution P p is much more interesting and complicated. 



As was shown in |£2| its behavior depends strongly on initial conditions. If the distribution P p 
was initially concentrated at < 0*, where 0* is some critical value of the field, then it moves 
towards small in the same way as P c , i.e. it cannot become stationary. On the other hand, 
if the initial value of the field is larger than 0*, the distribution moves towards larger and 
larger values of the field 0*, until it reaches the field P , at which the effective potential of the 
field becomes of the order of Planck density M p (we will assume M p = 1 hereafter), where the 
standard methods of quantum field theory in a curved classical space are no longer valid. 



Some further steps towards the solution of this problem were made by Nambu and Sasaki pi 
and Mijic p2| . Their papers contain many important results and insights. However, Mijic ]32] 
did not have a purpose to obtain a complete expression for the stationary distribution P p (0, £|%). 
The corresponding expressions were obtained for various types of potentials V(<p) in Unfor- 
tunately, according to ||31|| , the stationary distribution P p ( 
at ^> p , i.e. at V((p) ^> 1, where the methods used in |3T] are inapplicable. 



,t\x) is almost entirely concentrated 



We will continue this investigation by writing the system of stochastic equations for P p . These 
equations can be obtained from eqs. (0), ( p~8|) by adding the term 3HP P , which appears due 
to the growth of physical volume of the Universe by the factor 1 + 3iJ(0) dt during each time 
-g§, |28[: 



interval dt f|3l 



dPp 
dt 

OP, 
dt 



I ff 3/2 (X) d_ ( H^{ X ) dPp\ _ V\X) dP p 
2 2tt d X { 2tt d X ) SH( X ) d X 



d ^ 3/2 (0) d /# 3/2 (0) 
7u\ : ' 



V 2vr 



2tt 



Hp' 
3^(0) ^ 



+mx)p P 



+ 3#(0)P P 



(20) 



(21) 



To find solutions of these equations one must specify boundary conditions. The boundary condi- 
tions at the end of inflation follow from the conservation of the probability flux fl2F 



_IYL P i 

4 V p[ 







(22) 



X = 0e 



The last boundary condition is especially interesting, since it indicates that at least in the vicinity 
of x — 0e the probability distribution with respect to x looks as a square of the tunneling wave 
function. However, we have found solutions of (f20l), (|21~1) to be rather stable with respect to 
modification of these boundary conditions, whereas the conditions at the Planck boundary = P 
do play a very important role. 



Investigation of these conditions poses many difficult problems. First of all, our diffusion 
equations are based on the semiclassical approach to quantum gravity, which breaks down at 
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super-Planckian densities. Secondly, the shape of the effective potential may be strongly modified 
by quantum effects at ~ (fi p . Finally, our standard interpretation of the probability distribution 
P c and P p breaks down at V(4>) > 1, since at super-Planckian densities the notion of a classical 
scalar field in classical space-time does not make much sense. 



However, these problems by themselves suggest a possible answer. Inflation happens only in 
theories with very flat effective potentials. At the Planck density nothing can protect the effective 
potential from becoming steep. Hence, one may expect that inflation ceases to exist at <ft > 'Ppi 
which leads to the boundary condition 

P p (<j )p ,t\x)=P P (<P,t\x P )=0, (23) 

where V(<P P ) = V( Xp ) = 0(1). 



There is also another, much more general reason to expect that inflation kills itself as the 
potential energy density approaches the Planck density V ~ 1. Indeed, the amplitude of fluc- 
tuations of the scalar field generated during the typical time 5t = H~ x is given by H/2ir, and 
their typical wavelength at that time is 0{H^ 1 ). This means that the energy density associated 
with the gradients of these perturbations is of the order of H A ~ V 2 . Thus, in the domains with 
V > 1 the gradient energy density ~ V 2 becomes larger than the potential energy density V. 
This violates one of the basic assumptions necessary for inflation in domains with V > 1. One 
may expect also that large gradients of energy on the scale comparable to the scale of the horizon 
H^ 1 lead to creation of black holes rather than to the permanent self-reproduction of inflationary 
/i-regions. 

Of course, one may argue that all our considerations do not make sense at densities larger 
than the Planck density. When the energy density in any /i-region approaches the Planck density, 
it may no longer be described in terms of classical space-time and should be just thrown away 
from our consideration. In particular, its volume should not be considered as contributing to the 
total volume of the Universe. Thus, such domains should be neglected in our definition of P p . In 
this case the distribution P p for the field will stop moving towards higher values of and will 
approach a stationary regime when this distribution will be shifted towards <fi ~ <p p . 

What we are saying is even stronger. Even if one makes an attempt to consider the domains 
with V > 1 as a part of classical space-time, many parts of these domains drop out from the 
process of inflation. This means that the total volume of inflationary /i-regions cannot grow as 
fast as e 3Ht . 



We do not know which of these arguments, if any, will survive in the future theory of all 
fundamental interactions. However, all these arguments point out in the same direction: For a 
phenomenological description of stochastic processes in classical space-time one should impose 
boundary conditions of the type of fl23] ) which do not permit penetration of inflation deep into 
the realm of super-Planckian densities. As we will show in p8[, the exact form of these boundary 



conditions is not very important; most of them allow the same class of solutions as the boundary 
condition (|23|). Therefore in this paper we will use these boundary conditions, assuming for 



definiteness that they are imposed at V((p p ) = V(x P ) = 1- 
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One may try to obtain solutions of equations (E^), ( pTj ) in the form of the following series of 
biorthonormal system of eigenfunctions of the pair of adjoint linear operators (defined by the left 
hand sides of the equations below): 

oo 

P P (<Mx) = £e A ^ s (xR(0) • (24) 



s=l 



Indeed, this gives us a solution of eq. (^1|) if 

1 # 3/2 (x) d ( H^{ X ) d A V'(x) d 

\—^r^ x) ) ~ mxjdx^ six) + mx) ' Ux) = Xs Mx) ■ (25) 

and 

\h (^1 F^>)) + 1 H + • - * ■ <*> 

The orthonormality condition reads 

r<t>v 

^s(x)^j(x)dx = 8 s j ■ (27) 



In our case (with regular boundary conditions) one can easily show that the spectrum of Xj is 
discrete and bounded from above. Therefore the asymptotic solution for P p [<f>,t\x) (in the limit 
t — > oo) is given by 

P p (<Mx) =e Al > 1 ( X )7T 1 (0) • (l + 0(e-( Al - A2 )')) . (28) 
Here ipi(x) is the only positive eigenfunction of eq. (f25|), Ai is the corresponding (real) eigenvalue, 



and 7Ti(0) is the eigenfunction of the conjugate operator (BB) with the same eigenvalue Ai. Note, 
that Ai is the highest eigenvalue, Re (Ai — A2) > 0. This is the reason why the asymptotic 
equation (BSj) is valid at large t. We have found [133 that in realistic theories of inflation a typical 
time of relaxing to the asymptotic regime, At ~ (Ai — A2) _1 , is extremely small. Typically it is 
only about a few thousands Planck times, i.e. about 10 -40 sec. This means that the normalized 
distribution 

P p (<M| X )=e- Al 'P p (<Mx) (29) 
rapidly converges to the time-independent normalized distribution 

PMX) = P P (<P,t^ oo\ X ) = MX) ?ri(0) ■ (30) 

It is this stationary distribution that we were looking for. Because the growing factor e~ Al * is the 
same for all <fi (and x)-> one can use Pp instead of P p for calculation of all relative probabilities. 
In particular, P p ((j)\x) gives us the fraction of the volume of the Universe occupied by the field 
0, under the condition that the corresponding part of the Universe at some time in the past 
contained the field x- The remaining problem is to find the functions ipi{x) and 7Ti(0), and to 
check that all assumptions about the boundary conditions which we made on the way to eq. 
are actually satisfied. 
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We have solved this problem for chaotic inflation in a wide class of theories including the 
theories with polynomial and exponential effective potentials V{4>) and found the corresponding 
stationary distributions P5| . Here we will present some of our results for the theories |0 4 and 

V D e a<t> . 

Solution of equations (|25|) and (p6~D for ipi(x) an d 7T i(0) m the theory |0 4 shows that these 
functions are extremely small at ~ e and \ ~ Xe- They grow at large and \i then rapidly 
decrease, and vanish at = % = <ft p . With a decrease of A the solutions become more and 
more sharply peaked near the Planck boundary. (The functions ip\ and 7Ti for the exponential 
potential have a similar behavior, but they are less sharply peaked near <p p .) A detailed discussion 
of these solutions will be contained in [p8| . The eigenvalues Ai corresponding to different coupling 
constants A are given by the following table: 



A 


1 


KT 1 


KT 2 


1(T 3 


10" 4 


10~ 5 


10~ 6 


Ai 


2.813 


4.418 


5.543 


6.405 


7.057 


7.538 


7.885 



One can find also the second eigenvalue A2. For example, for A = 10~ 4 one gets A2 = 6.789. This 
means that for A = 10~ 4 the time of relaxation to the stationary distribution is At ~ (Ai — A2)~ 4 ~ 
4M" 1 ~ 10~ 42 seconds — a very short time indeed. 

Note that the parameter Ai shows the speed of exponential expansion of the volume filled by a 
given field 0. This speed does not depend on the field 0, and has the same order of magnitude as the 



speed of expansion at the Planck density. Indeed, Ai should be compared to 3H((f)) = 2y67rV A (0), 
which is equal to 2\^Qn at the Planck density. It can be shown |28j that in the limit A — > the 



eigenvalue Ai also becomes equal to 2\/67r ~ 8.681. The meaning of this result is very simple: in 
the limit A — > our solution becomes completely concentrated near the Planck boundary, and Ai 
becomes equal to 3if (0 P ). 

At first glance, independence of the speed of expansion of volume e Alt on the value of the field 
may seem counterintuitive. The meaning of this result is that the domain filled with the field 
gives the largest contribution to the growing volume of the Universe if it first diffuses towards 
the Planckian densities, spends there as long time as possible expanding with nearly Planckian 
speed, and then diffuses back to its original value 0. 

But what about the field (f) which is already at the Planck boundary? Why do the correspond- 
ing domains not grow exactly with the Planckian Hubble constant H((p p ) = 2^/6tt/3? It happens 
partially due to diffusion and slow rolling of the field towards smaller 0. However, the leading 
effect is the destructive diffusion towards the space-time foam with > p . One may visualize 
this process by painting white all domains with V(4>) < 1, and by painting black domains filled 
by space-time foam with V(<j)) > 1. Then each time H~ 1 (<j) p ) the volume of white domains with 
~ (f) p grows approximately e 3 times, but some 'black holes' appear in these domains, and, as a 



result, the total volume of white domains increases only e 3Al / 2 ^ 67r times. This suggests (by anal 
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ogy with ||35|| ) calling the factor df = 3Ai/2\/67r 'the fractal dimension of classical space-time', 
or 'the fractal dimension of the inflationary Universe'. (Note that df < 3 for A 7^ 0; for example, 
df = 2.6 for A = 10 -5 .) However, one should keep in mind that the fractal structure of the 
inflationary Universe in the chaotic inflation scenario in general is more complicated than in the 
new or old inflation and cannot be completely specified just by one fractal dimension |28 . 



The distribution P p (<f)\x) = ^i(x)vri(0) which we have obtained does not depend on time t. 
However, in general relativity one may use many different time parametrizations, and the same 
physics can be described differently in different 'times'. One of the most natural choices of time 
in the context of stochastic approach to inflation is the time r = In = / H((f>(x,t),t) dt 



26], |27|| . Here a (x, t) is a local value of the scale factor in the inflationary Universe. By using this 
time variable, we were able to obtain not only numerical solutions to the stochastic equations, 
but also simple asymptotic expressions describing these solutions. For example, for the theory 
both the eigenvalue Ai and the 'fractal dimension' df (which in this case refers both to the 
Planck boundary at <p p and to the end of inflation at e ) are given bydy = Ai~3 — 1.1 V% and 
the stationary distribution is 



3 11 

pmmx) ~ M~wm> W(x) + 0.4 " n) ' * """r* (3 ~ M)4? 

Note that the first factor coincides with the square of the tunneling wave function [[Jlj! This 
expression is valid in the whole interval from <p e to <p p and it correctly describes asymptotic 
behavior of P p (<j>, t\x) both at x ~ Xe and at x ~ Xv 

A similar investigation can be carried out for the theory V{<p) = V Q e a ^. The corresponding 
solution is 

Pp(Mx) - -p("s4)) ~ 1) ' d4 " : ) y " /2(0) ' (32) 

This expression gives a rather good approximation for P p (4>,t\x) for all 4> and x- 

The main result of our work is that under certain conditions the properties of our Universe 
can be described by a time-independent probability distribution, which we have found for theories 
with polynomial and exponential effective potentials. A lot of work still has to be done to verify 
this conclusion, see p8| . However, once this result is taken seriously, one should consider its 



interpretation and rather unusual implications. 

When making cosmological observations, we study our part of the Universe and find that in 
this part inflation ended about t e ~ 10 10 years ago. The standard assumption of the first models 
of inflation was that the total duration of the inflationary stage was At ~ 10~ 35 seconds. Thus 
one could come to an obvious conclusion that our part of the Universe was created in the Big 
Bang, at the time t e + At ~ 10 10 years ago. However, in our scenario the answer is quite different. 



16 



Let us consider an inflationary domain which gave rise to the process of self-reproduction 
of new inflationary domains. For illustrative purposes, one can visualize self-reproduction of 
inflationary domains as a branching process, which gives a qualitatively correct description of 
the actual physical process we consider. During this process, the first inflationary domain of 
initial radius ~ within the time H~ l ((j)) splits into e 3 ~ 20 independent inflationary 

domains of similar size. Each of them contains a slightly different field 0, modified both by 
classical motion down to the minimum of V((p) and by long-wavelength quantum fluctuations of 
amplitude ~ H/2tt. After the next time step which will be slightly different for each of 

these domains, they split again, etc. The whole process now looks like a branching tree growing 
from the first (root) domain. The radius of each branch is given by H -1 ; the total volume of 
all domains at any given time t corresponds to the 'cross-section' of all branches of the tree at 
that time, and is proportional to the number of branches. This volume rapidly grows, but when 
calculating it, one should take into account that those branches, in which the field becomes larger 
than (j) p , die and fall down from the tree, and each branch in which the field becomes smaller than 
4> e , ends on an apple (a part of the Universe where inflation ended and life became possible). 

One of our results is that even after we discard at each given moment the dead branches 
and the branches ended with apples, the total volume of live (inflationary) domains will continue 
growing exponentially, as e Al *. What is even more interesting, we have found that very soon 
the portion of branches with given properties (with given values of scalar fields, etc.) becomes 
time-independent. Thus, by observing any finite part of a tree at any given time t one cannot 
tell how old the tree is. 



To give a most dramatic representation of our conclusions, let us see where most of the apples 
grow. This can be done simply by integrating e Xlt from t = to t = T and taking the limit as 
T — > oo. The result obviously diverges at large T as A^f 1 e AlT , which means that most apples 
grow at an indefinitely large distance from the root. In other words, if we ask what is the total 
duration of inflation which produced a typical apple, the answer is that it is indefinitely large. 

This conclusion may seem very strange. Indeed, if one takes a typical point in the root domain, 
one can show that inflation at this point ends within a finite time At ~ 10" 35 seconds. This is 
a correct (though model-dependent) result which can be confirmed by stochastic methods, using 
the distribution P c ((j),At\x) P^] - How could it happen that the duration of inflation was any 
longer than 10 -35 seconds? 

The answer is related to the choice between P c and P p , or between roots and fruits. Typical 
points in the root domain drop out from the process of inflation within 10~ 35 seconds. The number 
of those points which drop out from inflation at a much later stage is exponentially suppressed, 
but they produce the main part of the total volume of the Universe. Note that the length of each 



particular branch continued back in time may well be finite |36| . However, there is no upper limit 
to the length of each branch, and, as we have seen, the longest branches produce almost all parts 
of the Universe with properties similar to the properties of the part where we live now. Since 
by local observations we can tell nothing about our distance in time from the root domain, our 
probabilistic arguments suggest that the root domain is, perhaps, indefinitely far away from us. 
Moreover, nothing in our part of the Universe depends on the distance from the root domain, 
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and, consequently, on the distance from the Big Bang. 

Thus, inflation solves many problems of the Big Bang theory. However, after making all kinds 
of improvements of this theory, we are now winding up with a model of a stationary Universe, 
in which the notion of the Big Bang loses its dominant position, being removed to the indefinite 
past. One should note, that for all practical purposes the end of inflation in our part of the 
Universe can be considered as a Big Bang (or, more precisely, as a Pretty Big Bang), after which 
our part of the Universe became hot and looked like an expanding ball of fire. However, the 
global structure of the inflationary universe cannot be described by this simple picture. 

From this point of view, inflationary theory cannot be considered as a part of the Big Bang 
theory or as its improved version. On the contrary, the Big Bang theory can be embedded into 
inflationary cosmology, and after this embedding the Big Bang cosmology provides an excellent 
description of the local structure of the inflationary Universe. 
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